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ABSTRACT 

We present precise constraints on the normalization of the power spectrum of mass 
fluctuations in the nearby universe, erg, as a function of the mean local matter density, 
f2 m . Using the observed local X-ray luminosity function of galaxy clusters from the 
extended BCS and REFLEX studies, a mass-luminosity relation determined from 
Chandra and ROSAT X-ray data and weak gravitational lensing observations, and 
the mass function predicted by the Hubble Volume simulations of Evrard et al. , we 
obtain er 8 = (0.508 ± 0.019) O:~(°- 253±a024 \ with £l m < 0.34 at 68 per cent confidence. 
The degeneracy between as and f2 m can be broken using Chandra measurements of the 
X-ray gas mass fractions in dynamically relaxed clusters. Using this information and 
including Gaussian priors on the mean baryon density of the universe and the Hubble 
constant, we obtain as = 0.695 ± 0.042 and f2 m = 0.287 ± 0.036, for an assumed flat 
ACDM cosmology (marginalized 68 per cent confidence limits). Our results are in good 
agreement with some recent studies based on the local X-ray temperature function of 
clusters, the redshift evolution of the X-ray luminosity and temperature functions of 
clusters, early results from the Sloan Digitized Sky Survey, the most recent results 
from studies of cosmic shear, and combined analyses of the 2dF galaxy redshift survey 
and cosmic microwave background anisotropies. 

Key words: cosmological parameters - X-rays: galaxies: clusters - gravitational 
lensing — large-scale structure of the universe — X-rays: galaxies: clusters 



1 INTRODUCTION 

The X-ray luminosity function of galaxy clusters in the 
nearby universe provides a powerful cosmological probe. The 
observed luminosity function, n(L), can be combined with 
a relation linking the observed X-ray luminosity and mass, 
and the mass function, n(M), predicted by simulations, to 
obtain tight constraints on the combination of cosmologi- 
cal parameters f2 m and erg, where S7 m is the mean matter 
density of the local universe and as is the root-mean-square 
(rms) variation of the density held smoothed by a top hat 
window function of size 8/i _1 Mpc. 

Observationally, the keys to such studies are precise de- 
terminations of the local X-ray luminosity function of clus- 
ters and the relation linking the observed X-ray luminosities 
and total masses. X-ray selection currently offers the best 
way to identify massive galaxy clusters, and the local X-ray 
luminosity function has now been precisely determined by 
the BCS (Ebeling et al. 1997; Ebeling et al. 2000) and RE- 
FLEX (Bohringer et al. 2002) studies. The flux-limited BCS 
and REFLEX samples, which are based on data horn the 
ROSAT All-Sky Survey (RASS; Triimper 1993), together 
include ~ 750 clusters and cover approximately two thirds 



of the sky. Recently, significant effort has also been invested 
into measuring the local temperature function of clusters 
{e.g. Markevitch 1998; Pierpaoli et al. 2001; Ikebe et al. 
2002), which offers a complementary method for determin- 
ing cosmological parameters. At present, however, the tem- 
perature function samples are significantly smaller than the 
combined BCS-plus- REFLEX luminosity function data set, 
and the analysis of the temperature function is complicated 
by the fact that these samples are selected according to X- 
ray flux as well as temperature, requiring the use of both 
mass-temperature and temperature-luminosity relations in 
the analysis. 

Recent years have also seen significant efforts di- 
rected towards a precise calibration of the 'virial' rela- 
tions linking the observed luminosities, temperatures and 
masses of galaxy clusters {e.g. Horner, Mushotzky & Scharf 
1999; Nevalainen, Markevitch & Forman 2000; Finoguenov, 
Reiprich & Bohringer 2001; Allen, Schmidt & Fabian 2001b; 
Sanderson et al. 2003). In particular, the launch of the Chan- 
dra X-ray Observatory has permitted the first precise mea- 
surements of the temperature and mass profiles of relaxed 
clusters from X-ray data. Using a combination of Chandra 
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and gravitational lensing data, Allen et al. (2001b) con- 
firmed that luminous, relaxed galaxy clusters follow the sim- 
ple scaling relations predicted by theory, but that the nor- 
malization of the observed mass-temperature relation mea- 
sured within r2soo (where the mean enclosed mass density 
is 2500 times the critical density of the universe at the red- 
shifts of the clusters) is approximately 40 per cent lower than 
predicted by standard adiabatic simulations. This highlights 
the likely importance of additional physics such as cooling 
and pre-heating in the intracluster gas (see also Pearce et 
al. 2000; Thomas et al. 2002; Voit et al. 2002; Muanwong et 
al. 2002). 

Theoretically, the primary requirement for cosmologi- 
cal studies using the observed luminosity and/or tempera- 
ture functions of clusters is a precise prediction of the mass 
function. This has now been achieved for flat ACDM (and 
rCDM) cosmologies using the Hubble Volume simulations 
of Jenkins et al. (2001) and Evrard et al. (2002). 

In this paper we present precise constraints on erg and 
fi m based on the observed local luminosity function of the 
most X-ray luminous clusters in the extended BCS (Ebel- 
ing et al. 2000) and REFLEX samples, and a new calibra- 
tion, using pointed Chandra and ROSAT X-ray observa- 
tions and weak gravitational lensing results, of the mass- 
luminosity relation linking the masses of clusters measured 
within r 2 oo to their total 0.1 — 2.4 keV ROSAT luminosi- 
ties. Having determined our combined constraint on as and 
Qm, we show that the degeneracy between these parameters 
can be broken using Chandra results on the X-ray gas mass 
fractions in the most dynamically relaxed clusters. Including 
Gaussian priors on the mean baryon density of the universe 
(Q b h 2 = 0.0205 ± 0.0018; O'Meara et al. 2001), the Hub- 
ble constant (h = 0.72 ± 0.08; Freedman et al. 2001), and 
a theoretical bias factor (b — 0.93 ± 0.05; Bialek, Evrard 
& Mohr 2001) relating the asymptotic baryon fraction in 
the most X-ray luminous clusters to the mean value for 
the universe as a whole, we obtain <T8 = 0.695 ± 0.042 and 
Sl m = 0.287 ± 0.036 (marginalized 68 per cent confidence 
limits for an assumed flat ACDM cosmology). We compare 
our results to other measurements based on the local num- 
ber density of clusters, evolution of the X-ray luminosity and 
temperature functions, the 2dF galaxy redshift survey, cos- 
mic microwave background anisotropies, and measurements 
of cosmic shear. 

Throughout this paper, a flat ACDM cosmology with 
S7a = 1 — fim is assumed. In order to facilitate a di- 
rect comparison with previous X-ray studies, results on 
the masses, X-ray luminosities and X-ray gas mass frac- 
tions of individual clusters are quoted for a Hubble pa- 
rameter h — ffo/100 kms -1 Mpc -1 = 0.5 or, equivalently, 
h 5 o = i?o/50kms _1 Mpc _1 = 1.0. 



moving matter density of the universe, O m (0) is the mean, 
present matter density in units of the critical density, and 
p c (0) = 3Hq/8tyG is the critical density at redshift zero. 

Using a spherical overdensity algorithm to measure the 
masses of clusters within radii r2oo , where the mean enclosed 
mass density is 200 times the critical density of the universe 
at the redshift of interest, Evrard et al. (2002) show that for 
a fiat ACDM cosmology, the mass fraction /(cr _1 ) can be 
written as 



/(a" 1 ) = AcxpHW^+Bn, 



(1) 



where, for fi m = 0.3 and z = 0, A = 0.22, B = 0.73 and e = 
3.86. Evrard et al. (2002) also provide simple interpolations 
for A,B and e for other values of S7 m and z. The differential 
number density of clusters with mass M at redshift z is 



dn(M,z) _ fia-^piz) 



dlna 



M 



(2) 



where p(z) = p (l + z) 3 is the mean mass density of the 
universe at redshift z. Following Viana & Liddle (1999), we 
write 



a(R, z) = a»(z) 



R 



8k- 1 Mpc 



-7(H) 



where 

7(fl) = (0.3r + 0.2) 



2.92 + lo gl 



R 



8ft" 1 Mpc 



(3) 



(4) 



and r is the shape parameter of the cold dark matter transfer 
function. Following Sugiyama (1995), we set 



exp 



-fi b (0) 



h fi b (0) 



0.5 Q m (0) 



(5) 



where To = 2.726K is the temperature of the cosmic mi- 
crowave background (Mather et al. 1994) and fib(0) = 
(0.0205 ± 0.0018)/i~ 2 is the mean, present-day baryon den- 
sity in units of the critical density (O'Meara et al. 2001). 
For this calculation, we set h = 0.72 and fix the primordial 
spectral index n = 1. The redshift-dependent quantity as(z) 
is related to its present value, <ts(0), by 



a 8 (z) = (T S (0) 



1 



g(n m (z)) _ 
ff(n m (o)) i + z 



(6) 



THEORY: THE PREDICTED LUMINOSITY 
FUNCTION OF GALAXY CLUSTERS 



where, for a flat ACDM universe 



1 209Slm(z) »m(z) 2 , .4/7' 

70 + 140 140 +ll ^ z > 



Jenkins et al. (2001) show that the predicted mass function g(Q m (z)) = —Q m (z) 
of galaxy clusters of mass M at redshift z can be written as 2 
a function of lnu _1 (M, z), where a 2 (M, z) is the variance 
of the linearly evolved density field smoothed by a spher- 
ical top-hat filter of comoving radius 7?, enclosing a mass (1 + z) 3 

M = 47ri? 3 p /3. Here p = Sl m (0)p c (0) is the mean co- ^m(z) = O m (0) - _ ^-p - (1 - z)3fim(0) • ( 8 ) 



and 
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e.g. Viana & Liddle (1996). Combining equations 2 and 3, 
we obtain 



An(M,z) _ jp(z) 



AM 



3M 2 



(9) 



In Section 3.2.4 we show that for clusters with X-ray lumi- 
nosities, L, exceeding 10 45 /ig 2 ergs -1 in the 0.1 — 2.4 keV 
ROSAT band, the mass M measured within r2oo can be re- 
lated to the X-ray luminosity by a power-law model of the 
form 



E(z)M = M 



E(z) 



(10) 



where the best fit values, uncertainties on Mo and a, and 
the intrinsic scatter in the relation, are determined from 
Chandra and ROSAT X-ray data and weak gravitational 
lensing results. The evolution parameter 



e(z) = (i + z)y/(i + zn m + n A /(i + . 



A ) 



(11) 



accounts for the expected cosmological scaling of the rela- 
tionship (e.g. Bryan & Norman 1998; Mathiesen & Evrard 
2001). Applying the chain rule to equations 9 and 10, we 
obtain 



Table 1. The observed, binned X-ray luminosity function of 



the most X-ray luminous (Lx 0. 



> 10 45 /igQ erg s 1 ) galaxy 



clusters from the extended BCS and REFLEX studies. Column 
2 gives the mean 0.1 — 2.4 keV X-ray luminosity for each bin, 
Lx, in 10 44 /i^" 2 ergs _1 . Error bars indicate the bin boundaries. 
Column 3 gives the number of clusters in each bin and column 
4 the comoving space density in /i| Mpc ~ 3 (10 44 ergs -1 A 
ACDM cosmology with f2 m = 0.3 and Q A = 0.7 is assumed. 



n c\us 



n(L) 
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17 


1.32 ±0.32 
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-11 


REFLEX 




20 


1.56 ±0.34 


X 


10 


-9 




16.271 4 ;" 


20 


4.44 ±0.98 


X 


io- 
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29.95l^ 2 7 


20 


1.77 ±0.39 


X 


10" 


-11 



provide the most powerful constraints on cosmological pa- 
rameters, relatively little information is lost by restricting 
ourselves to the largest systems. 

The binned X-ray luminosity functions for clusters 
with Lx, 0.1-2.4 > 10 45 h^Q erg s _1 from the northern ex- 
tended BCS (Ebeling et al. 2000) and southern REFLEX 
(Bohringer et al. 2002) studies are summarized in Table 1. 
The mean redshift of the BCS clusters in this luminosity 



range is z — 0.21. 



dn(L,z) 7p(z)a 
AL = 3M 



E(z) 



(12) 



The predicted differential luminosity function (i.e. the co- 
moving number density of clusters at redshift z with lumi- 
nosities in an interval dL around L) given by equation 12 
can be compared with the observed X-ray luminosity func- 
tion from the BCS and REFLEX studies to constrain the 
combination of cosmological parameters fl m (0) and t7s(0), 
hereafter referred to as fi m and as, respectively. 



3 OBSERVATIONS 

3.1 The observed X-ray luminosity function of the 
most luminous galaxy clusters in the RASS 

For this study, we concentrate on the local (z < 0.3) X- 
ray luminosity function and restrict ourselves to the most 
luminous clusters, with Lx > 10 45 /i^ 2 erg s _1 in the 0.1 — 
2.4 keV ROSAT band (for a flat ACDM cosmology with 
f2 m = 0.3 and SIa = 0.7). This selection is facilitated by the 
large sample size and well-determined selection functions of 
the BCS and REFLEX data sets. 

The restriction to high luminosities reduces systematic 
uncertainties by matching the luminosity range of the lu- 
minosity function data to the range over which the mass- 
luminosity relation has been calibrated (Section 3.2.4). At 
lower luminosities, the effects of pre-heating and cooling 
in the intracluster gas are expected to become important 
and cause the mass-luminosity and mass-temperature rela- 
tions to deviate from simple power-law forms (e.g. Cava- 
liere, Menci & Tozzi 1997). Since the most massive clusters 



3.2 The observed mass-luminosity relation for the 
most X-ray luminous galaxy clusters 

In determining the mass-luminosity relation we have used 
mass measurements obtained from Chandra X-ray observa- 
tions of dynamically relaxed clusters and weak gravitational 
lensing results drawn from the literature. X-ray luminosi- 
ties are determined from pointed ROSAT observations. In 
total, the sample used to define the mass-luminosity rela- 
tion includes 17 clusters with precise mass estimates and 
X-ray luminosities Lx, 0.1-2.4 > 10 45 h^ ergs -1 , spanning 
the redshift range 0.08 < z < 0.47. 



3.2.1 Chandra mass measurements 

The Advanced CCD Imaging Spectrometer (ACIS) on 
Chandra permits direct, simultaneous measurements of the 
X-ray gas temperature and density profiles and, via the hy- 
drostatic assumption, the total mass distributions in galaxy 
clusters. We have used Chandra to obtain precise mass mea- 
surements for a sample of ten of the most X-ray luminous, 
dynamically relaxed clusters identified from the RASS. The 
relaxed dynamical states of the clusters are demonstrated 
by their regular X-ray and optical morphologies, X-ray tem- 
perature maps and, in 6/10 cases, from the availability of 
consistent, independent mass measurements from gravita- 
tional lensing studies (see Section 4.4). 

* We make the comparison between the observed and predicted 
luminosity functions at z = 0.21, the mean redshift of the BCS 
clusters with Lx, 0.1-2. 4 > 10 45 erg s _1 . Shifting this redshift 
by ±0.05 does not significantly change the results. 
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Table 3. Summary of the Chandra mass measurements. Column 2 gives the evolution parameter, E(z), appropriate for each cluster. 
Columns 3 and 4 summarize the best-fitting NFW model parameters: the scale radius, r s (in h^ Mpc) and concentration parameter, c. 
Columns 5 and 6 give the virial radii, r200 (in Mpc) and masses, M200 (in 10 14 Mq). Error bars are 68 per cent confidence limits 
for a single interesting parameter, determined from the \ 2 grids. A flat ACDM cosmology with S7 m = 0.3 and S1a = 0.7 is assumed. 
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Table 2. Summary of the Chandra observations. The dates of the 
observations are given in column 3. Column 4 lists the exposure 
times in ks. 





z 


Date 


Exposure 


Abell 478 


0.088 


2001 Jan 27 


42.4 


PKS0745-191 


0.103 


2001 Jun 16 


17.9 


Abell 963 


0.206 


2000 Oct 11 


36.3 


Abell 2390 


0.230 


1999 Nov 07 


9.1 


Abell 2667 


0.233 


2001 Jun 19 


9.6 


Abell 1835 


0.252 


1999 Dec 12 


19.6 


Abell 611 


0.288 


2001 Nov 03 


36.1 


MS2137-2353 


0.313 


1999 Nov 18 


20.6 


RXJ1347-1145 


0.451 


2000 Mar 05/Apr 29 


18.9 


3C295 


0.461 


1999 Aug 30 


17.0 



The Chandra observations were made using the ACIS 
and the back-illuminated S3 detector between 1999 August 
30 and 2001 November 3. We have used the level-2 event 
lists provided by the standard Chandra pipeline processing. 
These lists were cleaned for periods of background flaring 
using the CIAO software package, resulting in the net expo- 
sure times summarized in Table 2. 

The data have been analysed using the methods de- 
scribed by Allen et al. (2001a, 2002b) and Schmidt et al. 
(2001; these papers present detailed mass analyses of Abell 
2390, RXJ1347-1145 and Abell 1835, respectively). In brief, 
concentric annular spectra were extracted from the cleaned 
event lists, centred on the peaks of the X-ray emission from 
the clusters. (For RXJ1347-1145, the data from the south- 
east quadrant of the cluster were excluded due to ongoing 
merger activity in that region; Allen et al. 2002b.) The spec- 
tra were analysed using XSPEC (version 11.0: Arnaud 1996) , 
the MEKAL plasma emission code (Kaastra & Mewe 1993; 
incorporating the Fe-L calculations of Liedhal, Osterheld & 
Goldstein 1995), and the photoelectric absorption models of 
Balucinska-Church & McCammon (1992; the absorbing col- 
umn density was included as a free parameter in the fits, alle- 
viating problems associated with uncertainties in the quan- 



Table 4. Summary of the mass results based on the Dahlc et al. 
(2002) weak lensing observations. Columns 2 and 3 summarize 
the redshifts and evolution parameters for the clusters. Column 4 
lists the virial masses, M200 ( m 10 14 Mq), determined from 
fits to the observed tangential shear profiles using NFW models 
with a fixed concentration parameter, c = 5. No correction for the 
effects of correlated substructure has been applied. Error bars are 
68 per cent confidence limits for a single interesting parameter. 
A ACDM cosmology with f2 m = 0.3 and Qa = 0.7 is assumed. 
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-6.2 


Abell 267 


0.230 


1.122 


15.5 


+4.4 
-3.9 


Abell 1576 


0.299 


1.165 


17.5 


+4.3 
-3.8 


Abell 1995 


0.320 


1.179 


20.1 


+4.6 
-4.2 


Abell 1351 


0.328 


1.184 


42.3 


+7.8 
-6.9 



turn efficiency of the detectors at low energies). Only data in 
the 0.5 — 7.0 keV energy range were used. The spectra for all 
annuli were modelled simultaneously in order to determine 
the deprojected X-ray gas temperature profiles, under the 
assumption of spherical symmetry. 

For the mass modelling, azimuthally averaged surface 
brightness profiles were constructed from background sub- 
tracted, flat-fielded images with a 0.984 x 0.984 arcsec 2 pixel 
scale (2x2 raw detector pixels) . When combined with the de- 
projected spectral temperature profiles, the surface bright- 
ness profiles can be used to determine the X-ray gas mass 
and total mass profiles in the clusters. ^ For this analysis, 

T The observed surface brightness profile and a particular pa- 
rameterized mass model are together used to predict the temper- 
ature profile of the X-ray gas. (We use the median temperature 
profile determined from 100 Monte-Carlo simulations. The outer- 
most pressure is fixed using an iterative technique which ensures 
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we have used an enhanced version of the image deprojec- 
tion code described by White, Jones & Forman (1997) with 
distances calculated using the code of Kayser, Helbig & 
Schramm (1997). 

We have parameterized the cluster mass (luminous plus 
dark matter) profiles using a Navarro, Frenk & White (1997; 
hereafter NFW) model with 



p(r) = 



p c (z)S c 



(r/r 3 ) (1 +r/r s ) 



(13) 



where p(r) is the mass density, p c {z) = 3H(z) 2 /8nG is the 
critical density for closure at redshift z, r B is the scale ra- 
dius, c is the concentration parameter (c = r20o/r s ) and 
S c = 200c 3 /3 [ln(l + c) - c/(l + c)]. The normalizations of 
the mass profiles can also be expressed in terms of an effec- 
tive velocity dispersion, a = \/50 r B cH(z) (with r B in units 
of Mpc and H(z) in kms" 1 Mpc" 1 ). 

The best-fit NFW parameter values and 68 per cent 
confidence limits are summarized in Table 3. This table 
also lists the 'virial' radii, r200, where the mean enclosed 
density is 200 times the critical density of the universe at 
the redshifts of the clusters, and the masses within these 
radii, M200. (The uncertainties on parameters are deter- 
mined directly from the \ 2 grids.) Note that the Chandra 
data only cover the central regions of the clusters out to 
radii 0.2 — 0.5 r2oo and thus some extrapolation of the mod- 
els, assuming that the NFW parameterization remains valid 
to r2oo, is required in calculating the virial masses. 



3.2.2 Weak lensing mass measurements 

In order to expand the sample of clusters used to construct 
the mass-luminosity relation, and allow certain tests of the 
fairness of this relation (Section 4.4) , we have also included 
data for clusters with precise mass measurements from wide 
field weak gravitational lensing studies. In particular, we 
have included data from the study of Dahle et al. (2002), 
who present aperture mass profiles for eight clusters with 
Lx, 0.1-2.4 > 10 45 h^ 2 ergs" 1 (Section 3.2.3) obtained from 
wide-field imaging with the University of Hawaii (UH) 2.24m 
telescope and UH8K camera. One of these clusters, Abell 
963, is also in our Chandra sample. Since the Chandra and 
lensing data for this cluster give consistent M200 results, 
but the Chandra data provide tighter constraints, we use 
the Chandra result in our default analysis. (The Dahle et 
al. 2002 mass measurement for Abell 963 is, however, used 
in our analysis of the weak lensing subsample, discussed in 
Section 4.4.) The weak lensing mass measurements made 



with the UH8K camera used a control annulus of 550 arcsec, 
which corresponds to 3.4 ftj^Mpc for a cluster at z = 0.3. 

Unlike X-ray mass measurements, which are based on 
the hydrostatic assumption, lensing mass measurements are 
independent of the dynamical state of the gravitating mat- 
ter. As a result, there are no restrictions in the Dahle et al. 
(2002) sample on the dynamical states of the clusters; sev- 
eral of the systems appear to be undergoing major merger 
events. The inclusion of the Dahle et al. (2002) clusters al- 
lows us both to examine the effects of dynamical activity 
on the mass-luminosity relation and assess whether the use 
of Chandra data for dynamically relaxed clusters is likely 
to bias our determination of cosmological parameters (see 
Section 4.4). 

We have used the aperture mass profiles presented by 
Dahle et al. (2002) to recover the mean tangential shear pro- 
files for the clusters and have fitted these with NFW mod- 
els. In general, the lensing data are unable to constrain both 
the concentration parameter and scale radius of the NFW 
models and so we have fixed c = 5 for this analysis, a typi- 
cal value for such massive clusters inferred from simulations 
(e.g. Navarro et al. 1997), and consistent with the Chandra 
results listed in Table 3. The masses of the clusters deter- 
mined from the Dahle et al. (2002) data are summarized in 
Table 4. 

In addition to the Dahle et al. (2002) UH8K data, ac- 
curate weak lensing mass measurements are also available 
for Abell 2390 (Squires et al. 1996) and RXJ1347-1145 (Fis- 
cher & Tyson 1997). In both cases the lensing mass mea- 
surements at r2oo are in good agreement with the Chandra 
results (Allen et al. 2001a, 2002b). Dahle et al. (2002) also 
present a weak lensing mass measurement for Abell 1835 us- 
ing a smaller camera, which is consistent with the Chandra 
result in Table 3. 

On the basis of numerical simulations, Metzler, White 
& Loken (2001) argue that large scale structure in the envi- 
ronments (r ~ 10 — 20/i" Mpc) of galaxy clusters are likely 
cause measurements of M200 from weak lensing to overesti- 
mate the true masses of clusters by, on average, ~ 30 per 
cent. Previous work by Cen (1997) and Reblinksy & Bartel- 
mann (1999) had argued for smaller effects, of the order of 
~ 10 per cent. Based on these studies, we have included a 
statistical correction of 20 per cent to the lensing masses in 
our determination of cosmological parameters i. e. we multi- 
ply the masses in Table 4 by 0.83. We note that the effects of 
more distant, uncorrelated structure along the lines of sight 
to the clusters are not expected to bias the lensing mass 
measurements (e.g. Metzler et al. 2001, Hoekstra 2002). 



a smooth pressure gradient in these regions.) The predicted tem- 
perature profile is rebinncd to the same binning as the spectral 
results and the \ 2 difference between the observed and predicted, 
deprojected temperature profiles is calculated. The parameters for 
the mass model are stepped through a regular grid of values in 
the r B -c or r B -cr planes to determine the best-fitting values and 68 
per cent confidence limits. (The best-fit models generally provide 
good descriptions of the data). The gas mass profile is determined 
to high precision at each grid point directly from the observed sur- 
face brightness profile and model temperature profile. Spherical 
symmetry and hydrostatic equilibrium are assumed throughout. 
Small departures from spherical symmetry are expected to have 
a negligible effect on the results e.g. Piffaretti et al. (2003). 



3.2.3 ROSAT luminosity measurements 

In constructing the mass-luminosity relation, we have used 
pointed ROSAT observations to determine the total, intrin- 
sic 0.1 — 2.4 keV luminosities of the clusters. This minimizes 
systematic uncertainties by matching the observing band 
and (as far as possible) detector technology to that used for 
the RASS observations. The agreement between the BCS 
and REFLEX fluxes, which are based on RASS data, and the 
fluxes determined from deep, pointed ROSAT observations 
is good (Ebeling et al. 1998; Bohringer et al. 2002). The de- 
tails of the pointed ROSAT observations are summarized in 
Table 5. The data were analysed using the XSELECT pack- 
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Table 5. Summary of the pointed ROSAT observations. Columns 2 and 3 list the date of observation and the detector used. Column 4 
gives the exposure times in ks. Column 5 lists the intrinsic 0.1 — 2.4 luminosities, Lo.i-2.4 (in 10 44 h^g ergs" 1 ). Error bars arc 68 per 
cent confidence limits. A ACDM cosmology with S7 m = 0.3 and Q A = 0.7 is assumed. 





Date 




Detector 


Exposure 


L 


1. 


L-2.4 


Abell 4 (8 


1991 


Aug 


31 




21.4 


22 


7 


±0 


2 


FKbO/4o-191 


1993 


Oct 


15 




1 a c: 
10. 


no 

28 


2 


±0 


6 


ADCll OzO 


1998 


Mar 


09 


Till T 

rirtl 


21 . I 


18 





±0 


8 


Abell 209 


1996 


Jul 


01 


HRI 


10.6 


15 


2 


± 1 





Abell 963 


1991 


Apr 


20 


HRI 


10.5 


13 


4 


± 1 





Abell 141 


1996 


Dec 


10 


HRI 


16.2 


12 


6 


±0 


7 


Abell 267 


1996 


Jan 


03 


HRI 


15.7 


11 


1 


±0 


9 


Abell 2390 


1993 


Nov 


13 


PSPC 


10.3 


31 


7 


±0 


5 


Abell 2667 


1994 


Dec 


14 


HRI 


21.3 


29 


2 


± 1 


1 


Abell 1835 


1993 


Jul 


03 


PSPC 


6.2 


38 


3 


±0 


9 


Abell 611 


1996 


Apr 


04 


HRI 


17.3 


11 


4 


± 1 


1 


Abell 1576 


1993 


Nov 


07 


PSPC 


16.3 


13 


4 


±0 


4 


MS2137-2353 


1993 


Nov 


07 


PSPC 


10.5 


19 


1 


±0 


9 


Abell 1995 


1995 


Nov 


13 


HRI 


16.5 


13 


7 


± 1 


4 


Abell 1351 


1995 


Apr 


29 


HRI 


31.9 


16 


5 


±3 





RXJ1347-1145 


1995 


Jan 


28 


HRI 


15.8 


81 


1 


±3 


1 


3C295 


1995 


Jun 


19 


HRI 


22.2 


12 


2 


± 1 


3 



"o 



X 



N 



in 



o 



X 
LO 



'88 



20 



50 



/E(z) [10 44 h -erg s 



"0.1-2.4 



- 1 ] 



Figure 1. The observed mass-luminosity relation. Chandra mass 
measurements for dynamically relaxed clusters are indicated by 
filled circles. Weak lensing mass measurements are indicated by 
open squares. A 20 per cent correction for the effects of corre- 
lated substructure has been applied to the lensing masses (Sec- 
tion 3.2.2). The four clusters with Chandra mass measurements 
and consistent weak lensing results are indicated by filled circles 
surrounded by open squares. The data for Abell 267 have been off- 
set slightly for display purposes. The best-fitting power law model 
from Section 3.2.4 is shown as the dashed curve. The two most 
significant outliers above (Abell 1351) and below (Abell 209) the 
best-fit curve appear to be undergoing major merger events (see 
Section 4.4). A ACDM cosmology with Q m = 0.3 and Q A = 0.7 
is assumed. 



age (version 2.0) and XSPEC (version 11.0). The emission- 
weighted temperatures and metallicities of the clusters were 
set to the values measured with Chandra. Where Chandra 
data were not available, the metallicity was set to 0.3 solar 
and the temperature was determined iteratively from the 
luminosity-temperature relation of Allen & Fabian (1998). 
Note, however, that the precise settings of the temperatures 
and metallicities have little effect on the measured 0.1 — 2.4 
keV luminosities for such hot, massive clusters. The absorb- 
ing column densities were set to the Galactic values deter- 
mined by Dickey & Lockman (1990). 

The intrinsic 0.1 — 2.4 keV luminosities for a ACDM 
cosmology with fi m = 0.3, = 0.7 and h = 0.5 (the same 
cosmology assumed in the BCS and REFLEX luminosity 
functions in Table 1) are listed in Table 5. 

3.2.4 The observed mass-lurninosity relation 

The mass-luminosity relation for the 17 X-ray luminous clus- 
ters in our sample, measured within radii r2oo corresponding 
to a density contrast A = 200 with respect to the critical 
density of the universe at the redshifts of the clusters, is 
shown in Fig. 1. Those clusters with mass measurements 
from Chandra X-ray data are indicated by filled circles. 
Clusters with weak lensing mass measurements are marked 
with open squares. The four clusters with Chandra mass 
measurements and consistent weak lensing results are indi- 
cated by filled circles surrounded by open squares. A flat 
ACDM cosmology with fi m = 0.3 and J1a = 0.7 is assumed. 
We have fitted the data using a model of the form 



logio 



E(z) M200 



= a login 



+ logi 



E(z) lO 44 /^ 2 ergs- 
Mo 



(14) 
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a 

Figure 2. The 68.3 and 95.4 per cent confidence contours on the 
intercept and slope of the model fitted to the mass-luminosity 
data. The plus sign marks the best-fit position. 

Using the BCES(F|X) estimator of Akritas & Bershady 
(1996), which accounts for errors in both axes and the 
presence of possible intrinsic scatter, we obtain best-fitting 
values and 68 per cent confidence limits from 10 6 boot- 
strap simulations of logio(M //i 50 1 M Q ) = 14.29^Q 23 and 
a — 0.761q 13 \ The distribution of Mo and a values are 
shown in Fig. 2. The observed slope is in good agreement 
with the expected slope of the mass-bolometric luminosity 
relation of a — 0.75, from models of simple gravitational 
collapse. 

The rms scatter of the observed log [E(z)M2oo] values 
about the best fitting curve is 0.22 for the full sample of 17 
clusters, 0.15 for the 10 dynamically relaxed clusters stud- 
ied with Chandra, and 0.29 for the 8 clusters with weak 
lensing measurements from Dahle et al. (2002). The scatter 
in the full sample is similar to that measured by Reiprich 
& Bohringer (2001) using ASCA and ROSAT X-ray data. 
Note, however, that our best-fitting curve implies a slightly 
higher X-ray luminosity for a given mass. 

4 DETERMINATION OF COSMOLOGICAL 
PARAMETERS 

4.1 Monte Carlo method 

Using the theoretical prescription for the mass function of 
galaxy clusters described in Section 2, the observed BCS and 
REFLEX X-ray luminosity functions summarized in Table 
1, and the mass-luminosity data discussed in Section 3.2.4, 
we can determine as as a function of f2 m . 

In determining our results on cosmological parameters, 
we have used a Monte Carlo approach. For each iteration of 
the code, we construct a random bootstrap sample of the 



1.4 



b 




0.4 - 



0.2 0.4 0.6 0.8 1 



Figure 3. The 68.3 and 95.4 per cent confidence contours on <rg 
and f2 m obtained from the BCS+REFLEX luminosity function 
data (Tabic 1), theoretical mass function, and observed mass- 
luminosity relation (Section 3.2.4). A flat ACDM cosmology is 
assumed. 



M200 and Lo. 1-2.4 values listed in Tables 3 — 5 and examine 
a grid of (96 x 131) fi m and as values, covering the plane 
0.05 < ri m < 1.0 and 0.20 < a 8 < 1.50. For each value of 
Qm, we scale the mass and luminosity values appropriately 
and determine the best-fitting mass-luminosity relation. For 
each f2 m , a 8 parameter pair, we construct a model X-ray 
luminosity function, including the effects of random scatter 
in the mass- luminosity relation, which we characterize us- 
ing a log-normal distribution. The model X-ray luminosity 
function is then compared with the observed BCS and RE- 
FLEX data, also scaled to the appropriate cosmology, and 
the \ 2 difference between the two is calculated. In this way, 
the best fitting f2 m , as pair for the grid is determined, which 
provides us with a single sample result. The whole process 
is repeated for 10 6 iterations to produce the final results, 
discussed below. 



4.2 Results on as as a function of fi m 

Fig. 3 shows the 68.3 and 95.4 per cent confidence con- 
tours in the as — f2 m plane from one million iterations of 
the Monte Carlo code. The results exhibit the well known 
degeneracy between as and Q m , which can be approxi- 
mated (for 0.1 < fi m < 0.4) by the simple fitting formula 
as = (0.510 ± 0.019) q~(°- 2B3±0 - 024 ) _ Q ur ana lysis favours 
low values for f2 m : marginalizing over as, we find fi m < 0.34 
at 68 per cent confidence. 

Note that when neglecting the effects of scatter and 
uncertainties in the normalization and slope of the mass- 
luminosity relation, we obtain the best fit (\ 2 = 5.1 for 4 
degrees of freedom) for f2 m = 0.23 and as = 0.74. 
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4.3 Breaking the erg — f2 m degeneracy using the 
Chandra f gas (z) data 

The usual approach, discussed in the literature, to break the 
degeneracy between erg and il m (Fig. 3) is to use the red- 
shift evolution of the luminosity and/or temperature func- 
tion of clusters, which depends strongly on the mean mass 
density of the universe (see references in Section 5.2). How- 
ever, this approach is both observationally challenging in 
terms of identifying complete, high-redshift cluster samples, 
and prone to systematic uncertainties due to potentially 
increased levels of dynamical activity and contaminating 
AGN emission at high redshifts, which introduce additional 
scatter into the mass-luminosity and mass-temperature re- 
lations. 

Fortunately, the Chandra data offer a powerful alter- 
native method to break the degeneracy between erg and fi m 
using the observed X-ray gas mass fractions, / gas , in the clus- 
ters and their apparent redshift dependence (e.g. White & 
Frenk 1991; White et al. 1993; Sasaki 1996; Pen 1997; Ettori 
& Fabian 1999; Allen et al. 2002a; Erdogdu, Ettori & Lahav 
2002; Ettori, Tozzi & Rosati 2003). The matter content of 
rich clusters of galaxies is thought to provide an almost fair 
sample of the matter content of the universe [e.g. White et 
al. 1993). The observed ratio of baryonic to total mass in 
clusters is therefore expected to closely match the ratio of 
the cosmological parameters S7b/^m, where fib is the mean 
baryon density of the universe. The apparent redshift de- 
pendence of the / gas measurements arises from the fact that 
the measured / gas values depend upon the assumed angular 
diameter distances to the sources as / gas oc D A 5 . Thus, al- 
though we expect the measured / gas values to be invariant 
with redshift, this will only appear to be the case when the 
assumed cosmology matches the true, underlying cosmology. 

The observed / gas profiles for the ten relaxed clusters 
studied with Chandra, for an assumed h = 0.5 standard 
cold dark matter (SCDM) cosmology with fi m = 1.0 and 
SIa = 0.0, are shown in Fig. 4. (The six clusters previously 
studied by Allen et al. 2002a are shown in a lighter shad- 
ing.) With the possible exception of Abell 963, we see that 
the / gas profiles appear to have converged, or be close to 
converging, within r2soo- We note that the data for the two 
nearest clusters, Abell 478 and PKS0745-191, do not extend 
to r25oo- However, their / gas profiles appear to be close to 
converging within r ~ 0.6r25oo- 

Fig. 5 shows the / gas measurements at r25oo (or at the 
outermost radii studied for Abell 478 and PKS0745-191) as 
a function of redshift, for the nine clusters with convergent 
/ g a S profiles. Following Allen et al. (2002a), we have fitted 
these data with the model 



rmod/ s 
J g as 



(l+0.19Vh) fin 



2j9» =1 ' nA=0 (*) 



0.5 D? n 



, n A =i-n„ 



(15) 



and determined the best-fitting value of f2 m , for an assumed 
fiat ACDM cosmology. The parameter b is a bias factor 
that is motivated by gasdynamical simulations, which sug- 
gest that the baryon fraction in clusters is slightly depressed 
with respect to the universe as a whole (e.g. Cen & Ostriker 
1994; Eke, Navarro & Frenk 1998a; Frenk et al. 1999; Bialek 
et al. 2001). We include a Gaussian prior on the bias factor, 
b = 0.93 ± 0.05, a value appropriate for hot (kT > 5 keV), 




2500 



Figure 4. The X-ray gas mass fraction profiles measured with 
Chandra, with the radial axes scaled in units of r2soo . The results 
for the six clusters previously discussed by Allen et al. (2002a) 
are shown in lighter shading. The new results for Abell 478, 611, 
963 and 2667 are shown as dark circles. Note that / gas (r) is an 
integrated quantity and so the error bars on neighbouring points 
in a profile are correlated. An SCDM cosmology with h = 0.5 is 
assumed. 



massive clusters in the redshift range < z < 0.5 from the 
simulations of Bialek et al. (2001). We also include Gaus- 
sian priors on the Hubble constant, h = 0.72 ± 0.08, the 
final result from the Hubble Key Project reported by Freed- 
man et al. (2001), and Q h h 2 = 0.0205 ± 0.0018 (O'Meara 
et al. 2001), from cosmic nucleosynthesis calculations con- 
strained by the observed abundances of light elements at 
high redshifts. The constraints on fl m determined from this 
analysis are shown as the dark, solid curve in Fig. 6. We ob- 
tain fi m = 0.291±o;°36 at 68 per cent confidence. Also shown 
(dashed curve) are the results obtained when fixing the bias 



parameter b = 1.0, for which fi m = O.SMto^ . 

We can now combine (by multiplying the relevant prob- 
ability densities) the constraints on f2 m from Fig. 6 with the 
joint constraints on erg and Q m , shown in Fig. 3, to obtain 
our final results on erg and Q m . These are shown in Fig. 7. 
Using the Gaussian prior on the bias parameter, we obtain 
a 8 = 0.695 ± 0.042 and fi m = 0.287 ± 0.036 (marginalized 
68 per cent confidence limits). With b = 1.0 fixed, we obtain 
erg = 0.683 ± 0.041 and fi m = 0.309 ± 0.035. 



4.4 Systematic uncertainties and the effects of 
merger events 

An important aspect of the present work is the reduced 
level of systematic uncertainty with respect to most pre- 
vious studies based on the local abundance of galaxy clus- 
ters. In the first case, the independent BCS and REFLEX 
studies have well determined selection functions (Ebeling et 
al. 1998, 2000; Bohringer et al. 2002) and provide precise, 
consistent results on the local X-ray luminosity function. 
With the large size of the combined BCS-plus-REFLEX data 
set (which covers two thirds of the sky) we have been able 
to limit our analysis to the most luminous clusters, with 
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Figure 5. The apparent redshift variation of the X-ray gas mass 
fraction measured at T2500 (with root-mean-square ltr errors) for 
the nine clusters with convergent / gas profiles in Fig. 4 (see text). 
The results for the six clusters previously discussed by Allen et 
al. (2002a) are shown in lighter shading. The solid curve shows 
the predicted /gas(z) behaviour for a flat ACDM cosmology with 
Q m = 0.291. 



O 



O 



o 




0.1 0.2 0.3 0.4 0.5 



Figure 6. The constraints on n m from the Chandra / gas (^) data 
in Fig. 5. A flat ACDM cosmology and Gaussian priors of h = 
0.72 ± 0.08 and tt h h 2 = 0.0205 ± 0.0018 are assumed. The solid 
curve shows the results obtained using a Gaussian prior on the 
bias factor, b = 0.93±0.05. The dashed curves show the results for 
b = 1 (fixed). The 1, 2 and 3 sigma confidence limits are marked 
as dotted lines. 
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Figure 7. The 68.3, 95.4 and 99.7 per cent confidence contours 
on <rg and n m from the combined analysis of the BCS+REFLEX 
luminosity function and Chandra f gaB (z) data. A flat ACDM cos- 
mology is assumed. 

L x ,o.i-2A > lO^/i^ergs" 1 (for an fi m = 0.3, = 0.7, 
cosmology) , for which the mass- luminosity relation has been 
calibrated using Chandra and ROSAT X-ray data and weak 
lensing observations (Section 3.2.4). The complicating ef- 
fects of cooling and pre- heating are minimized for such mas- 
sive clusters and a simple power-law model provides a rea- 
sonable description of the mass- luminosity relation (Fig. 1), 
albeit with some residual intrinsic scatter that is in part 
related to the different dynamical histories of the clusters 
and residual systematic effects in the mass measurements 
(Sections 3.2.2, 3.2.4; see also below). Repeating the analy- 
sis using only the BCS luminosity function data, or only the 
REFLEX data, or using a higher luminosity cut in the lu- 
minosity function (e.g. Lx,o.i-2A > 1.3 X 10 45 ergs" 1 ), 
leads to consistent results (although with larger statistical 
uncertainties). Consistent results on cosmological parame- 
ters are also obtained if we replace the Evrard et al. (2002) 
critical spherical overdensity mass function (Section 2) with 
the Jenkins et al. (2001) mean spherical overdensity (S0324) 
prescription, scaling the cluster masses accordingly. 

The dominant statistical and systematic uncertainties 
in the analysis are associated with the mass-luminosity rela- 
tion. However, care has been taken to minimize the sys- 
tematic uncertainties. In particular, we have limited our 
study of the mass-luminosity relation to clusters with pre- 
cise mass measurements from Chandra or wide-field grav- 
itational lensing studies. The clusters studied with Chan- 
dra have been selected from the RASS as the most lumi- 
nous, dynamically relaxed (in terms of their X-ray and op- 
tical morphologies) clusters known. The relaxed natures of 
these systems means that they are the clusters for which X- 
ray mass measurements are most reliable. For Abell 963, 
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1835, 2390 and RXJ1347-1145, independent confirmation 
of the Chandra mass measurements at r2oo is available 
from weak gravitational lensing studies (Dahle et al. 2002; 
Squires et al. 1996; Fischer & Tyson 1997; Allen et al. 
2001a, 2002b). A programme of weak lensing measurements 
for the other clusters in our Chandra sample is underway 
(Gray et al. , in preparation). In addition, consistent strong 
lensing mass measurements are available for Abell 1835, 
2390, RXJ1347.5-1145, MS2137.3-2353 and PKS0745-191 
(Schmidt et al. 2001, Allen et al. 2001a, 2002b, Schmidt et 
al. , in preparation) . The presence of significant non-thermal 
pressure support on large spatial scales in these clusters can 
therefore be excluded. 

As well as increasing the size of our calibration sam- 
ple, the inclusion of the Dahle et al. (2002) weak lensing 
data is important in that it allows us to examine the dis- 
tribution of less dynamically relaxed systems in the mass- 
luminosity plane. Whereas X-ray data can only be used to 
obtain precise mass measurements for dynamically relaxed 
clusters (the X-ray mass measurements are based on the as- 
sumption of hydrostatic equilibrium in the X-ray gas), lens- 
ing mass measurements are essentially independent of the 
dynamical state of the gravitating matter and so can be 
extended to systems undergoing merger events. The Dahle 
et al. (2002) sample includes clusters with a range of dy- 
namical states. Inspection of ROSAT and Chandra images 
of these systems shows that many exhibit significant dy- 
namical activity, with Abell 141, 209, 520, 1351 and 1576 
undergoing major merger events. (Similar conclusions are 
drawn by Dahle et al. 2002 from their optical data). The 
BCS and REFLEX samples, from which the X-ray luminos- 
ity functions have been constructed (Section 3.1), include 
all clusters above the respective X-ray flux limits, indepen- 
dent of their dynamical states. It is therefore important that 
the mass-luminosity relation provides a fair translation be- 
tween luminosity and mass for all clusters included in the 
samples. Although the 17 clusters in our mass-luminosity 
relation do not represent a complete subsample of the BCS 
and REFLEX data sets, inspection of Fig. 1 suggests that 
the Dahle et al. (2002) clusters exhibit a similar mean mass 
per unit X-ray luminosity to the relaxed systems studied 
with Chandra, albeit with two significant outliers from the 
best-fit curve (see below) and, therefore, that the applica- 
tion of the mass-luminosity relation determined here to the 
BCS and REFLEX samples is reasonable. 

We have carried out a more rigorous test of the va- 
lidity of the mass-luminosity relation by repeating the de- 
termination of cosmological parameters using other mass- 
luminosity relations constructed from subsamples of the cal- 
ibration data. In the first case, a mass- luminosity relation 
based on the Chandra mass measurements (10 clusters) 
was examined. In the second case, only the weak lensing 
measurements were used (8 clusters). * These results were 
then compared with those obtained from the full sample 
(17 clusters, excluding the lensing result for Abell 963). 
The normalized, marginalized probability distributions for 

t For the determination of cosmological parameters with the 
weak lensing subsample, only the first two bins of the BCS and 
REFLEX luminosity functions (Table 1) were fitted, as the Dahle 
et al. (2002) clusters do not populate the highest luminosity bin. 




0.4 0.6 0.8 1 1.2 

Figure 8. The normalized, marginalized probability distribu- 
tions for ag obtained from the Monte Carlo analysis using mass- 
luminosity relations constructed from the whole sample of 17 clus- 
ters (solid line), the subsample of 10 clusters with Chandra mass 
measurements (dotted line), and the subsample of 8 clusters with 
wide-field weak lensing mass measurements (grey line). See text 
for details. 



as are shown in Fig. 8. For the Chandra subsample, we ob- 
tain as = 0.69 ± 0.04. For the lensing subsample, we find 
as = 0.71^0 12 (68 per cent confidence limits). Both re- 
sults are consistent with those obtained for the whole sam- 
ple, as = 0.70 ± 0.04, although for the lensing subsample 
the increased scatter in the mass-luminosity relation leads 
to larger error bars. It appears, then, that the inclusion of 
Chandra mass measurements for dynamically relaxed clus- 
ters does not bias the best-fit result on as significantly, but 
does reduce the formal statistical uncertainties. Assessing 
whether this leads us to underestimate the true uncertainties 
in as will require more data and the ability to disentangle 
whether the scatter in the lensing results arises primarily 
from the effects of correlated substructure (or similar sys- 
tematic problems in the lensing analysis), or the effects of 
dynamical activity on the X-ray gas (see below). 

Ricker & Sarazin (2001), Ritchie & Thomas (2002) and 
Randall, Sarazin & Ricker (2002) present simulations of the 
effects of mergers on the X-ray properties of clusters for a 
variety of mass ratios, impact parameters and central gas 
densities. These authors suggest that major mergers can 
lead to significant short-term increases in the luminosities 
of clusters during the periods of closest approach, which are 
then followed by dips in luminosity as the merging dark 
matter cores move apart, before the cluster returns to equi- 
librium. From Fig. 1 we see that Abell 209 appears to have 
an unusually high X-ray luminosity /mass ratio, which may 
have been boosted by the ongoing merger activity in this 
cluster. (We note, however, that MS2137. 3-2353, which is a 
highly relaxed cluster with no obvious merger activity and 
a sharp central density peak, also appears to have a high 
luminosity /mass ratio and lies below the best-fitting curve. 
In this case, the offset may indicate a relatively early for- 



© 0000 RAS, MNRAS 000, 000-000 



Cosmological constraints from the local X-ray luminosity function of galaxy clusters 11 



mation epoch for the cluster.) In contrast Abell 1351, which 
is also undergoing a major merger event, has an unusually 
low X-ray luminosity /mass ratio and lies above the best-fit 
curve in Fig. 1. Detailed simulations and further observa- 
tions are required to improve our understanding of the ef- 
fects of mergers on the X-ray properties of clusters. However, 
the indications from the present study are that the effects of 
mergers on the mean mass-luminosity relation for the most 
luminous clusters are relatively small. 

Chandra measurements of the X-ray gas mass fraction, 
/g as , in dynamically relaxed clusters provide one of the most 
simple and robust methods by which to measure fi m (see 
e.g. Allen et al. 2002a and references therein). The largest 
systematic uncertainties in the analysis lie in how well the 
measured baryonic mass fractions in the clusters approxi- 
mate the universal mean. Lowering the bias factor, b, by 
~ 10 per cent would cause the best-fitting value of fi m to 
rise by a similar amount. Similarly, although the / gas profiles 
in Fig. 4 appear to have converged, or be close to converging, 
at the outer measurement radii, any rise in the / gas values 
beyond these points would cause a corresponding reduction 
in fi m . Any change in the best-fitting value of fi m would 
then lead to a change in as, as shown in Fig. 3. 

Finally, we note that our Monte Carlo analysis (Section 
4.1) takes account of the intrinsic scatter and uncertainties 
in the normalization and slope of the mass-luminosity rela- 
tion. (The scatter is modelled as a log-normal distribution 
about the best-fitting curve, which provides a good descrip- 
tion of the measured offsets.) Neglecting such effects would 
cause to the best-fitting value of as to rise slightly, and sig- 
nificantly reduce the formal statistical uncertainties on the 
measured parameters. 

5 COMPARISON WITH OTHER RESULTS 

5.1 Other local cluster abundance studies 

Fig. 9 shows a comparison of the results on erg as a func- 
tion of fi m from the present study (thick, solid curves: as 
in Fig. 3), with the findings from five other recent studies 
based on the observed local number density of galaxy clus- 
ters. The dot-dashed curve in Fig. 9 shows the result of Sel- 
jak (2002), a 8 = (0.44 ± 0.04) fi m °' 44 (r/0.2) ' 08 , using the 
observed local X-ray temperature function of rich clusters 
(Pierpaoli, Scott & White 2002) and a mass-temperature 
relation determined from ASCA and ROSAT observations 
(Finogucnov et al. 2001). The result of Seljak is consis- 
tent with ours for fi m ~ 0.3, although the present study 
leads to tighter constraints. The result of Pierpaoli et al. 
(2001), a 8 = (0.50 ± 0.04) fi m 0(i0 (dotted curve), based 
on the same local temperature function data but normal- 
ized by a theoretical mass-temperature relation, lies sig- 
nificantly above ours. The result of Reiprich & Bohringer 
(2001), a s = 0.43 fi m ' 38 (short dashes), based on ASCA and 
ROSAT observations of RASS selected clusters, is in good 
agreement with the present study for fi m ~ 0.3. The result 
of Viana, Nichol & Liddle (2002), cr 8 = 0.38 fi m °' 48+0 ' 27 n ™ 
(long dashes), using the local REFLEX X-ray luminosity 
function and a mass-luminosity relation determined from 
ROSAT X-ray observations and a stacked weak lensing anal- 
ysis of relatively low-mass clusters identified in Sloan Digi- 
tized Sky Survey (SDSS) commissioning data, is consistent 



with the present study for fi m ~ 0.3. The result of Bah- 
call et al. (2002), a» = 0.35 fi m ' 60 , obtained by combining 
the number density of optically-selected, relatively low-mass 
clusters observed in SDSS commissioning data with a mass- 
optical richness correlation, is in good agreement with this 
work for fi m ~ 0.3. Finally, our results are in good agree- 
ment with the recent findings of Schuecker et al. (2003; not 
plotted), a» = 0.711±o.S|?, fim = 0.34ltg ^g, from a com- 
bined analysis of the X-ray luminosity function and large- 
scale clustering in the REFLEX sample. 



5.2 Evolution in the X-ray luminosity and 
temperature functions 

Borgani et al. (2001) present constraints on as and fi m de- 
termined from an analysis of evolution in the X-ray luminos- 
ity function of clusters in the ROSAT Deep Cluster Survey, 
which spans the redshift range z < 1.3. Their results of 
a 8 = 0.66+Q 05 and fi m = 0.35l°:io for a flat ACDM cos- 
mology are consistent with those reported here. 

Donahue & Voit (1999) report results from an anal- 
ysis of evolution in the temperature function of clusters 
within z < 0.8, using the low-redshift cluster sample of 
Markevitch (1998) and a high redshift sample identified 
from the Einstein Observatory Extended Medium Sensi- 
tivity Survey (EMSS). Their results, for an assumed flat 
ACDM cosmology, of a 8 = 0.73±g;g| and fi m = 0.27 ± 0.10 
are in good agreement with ours. Eke et al. (1998b) ob- 
tain as = 0.75 ± 0.15 and fi m = 0.36 ± 0.25 (flat ACDM) 
from an analysis of the temperature function within z ~ 0.4, 
which is consistent with the present work. Our results are 
marginally consistent with the values of as = 0.72 ±0.10 and 
fi m = 0.49 ± 0.12 reported by Henry (2000) from an analy- 
sis combining the local temperature function of clusters with 
the properties of EMSS clusters within z < 0.6. 



5.3 Cosmic shear measurements 

Fig. 10 shows a comparison of the results on as as a func- 
tion of fi m determined from the present study (dark, solid 
curves: as in Fig. 3) together with the findings from seven 
recent studies based on measurements of weak gravitational 
lensing due to large scale structure (cosmic shear). The up- 
per dotted curve in Fig. 10 shows the result of Van Waerbeke 
et al. (2002), a s = (0.57±0.04) n(°-24±o.i8)fi m -o.49_ The up _ 

per short-dashed curve shows the result of Refregier, Rhodes 
& Groth (2002), as = (0.55 ±0.08) fi m ' 44 . The long-dashed 
curve shows the result of Bacon et al. (2002), as = (0.43 ± 
0.06) fi m ' 68 . The dot-dashed curve shows the result of Hock- 
stra, Yee & Gladders (2002), a s = (0.45 ± 0.05) fi m ' 55 . The 
results on as from these four studies lie 20 — 35 per cent 
above the present work for fi m ~ 0.3. (Note that in order 
to obtain agreement with these results, the mean mass per 
unit X-ray luminosity for the clusters included in the present 
study would need to be raised by a factor ~ 2.5, well beyond 
the systematic uncertainties, which are < 20 per cent.) 

The lower three curves in Fig. 10 show the results from 
the three most recent cosmic shear studies, which have ap- 
peared on preprint servers after this paper was submitted. 
The lower dotted curve shows the result of Hamana et al. 
(2002), as = (0.41 ± 0.09) fi m ' 43 . The lower dashed curve 



© 0000 RAS, MNRAS 000, 000-000 



12 S. W. Allen et al. 




Figure 9. The 68.3 and 95.4 per cent confidence contours on as as 
a function of Q m determined from the present study (thick, solid 
curves; as in Fig. 3) together with the best-fit results of Seljak 
(2002: dot-dashed curve), Pierpaoli et al. (2001: dotted curve), 
Reiprich & Bohringer (2001: short-dashed curve), Viana et al. 
(2002: long-dashed curve) and Bahcall et al. (2002: thin, solid 
curve). See text for details. A flat ACDM cosmology is assumed. 

shows the result of Brown et al. (2002), ag = (0.41 ± 
0.05) ft" 50 . Finally, the solid curve shows the result of 
Jarvis et al. (2002), ct 8 = (0.36 ± 0.04) fC, 0,57 . Importantly, 
the Brown et al. (2002) and Jarvis et al. (2002) studies in- 
clude improved estimates of the redshifts of the lensed, back- 
ground sources. (A flat ACDM cosmology is assumed in all 
cases.) The results from the three most recent cosmic shear 
studies are in good agreement with the present work. 



5.4 Cosmic microwave background anisotropies 
and the 2dF galaxy redshift survey 

Fig. 11 shows a comparison of the 68.3 and 95.4 per cent 
confidence constraints on as and S7 m obtained from the 
present study (inner, shaded contours), including the Chan- 
dra f sas (z) data and Gaussian priors on the bias parameter 
(b = 0.93 ±0.05), Hubble constant (h = 0.72 ±0.08) and fl h 
(Q b h 2 = 0.0205 ± 0.0018), with the results from analyses of 
cosmic microwave background (CMB) anisotropies and 2dF 
galaxy redshift survey data. For this comparison, we have 
used the published Markov Chain Monte Carlo samples of 
Lewis & Bridle (2002) . The CMB data consist of a combina- 
tion of COBE (Bennett et al. 1996), Boomerang (Netterfield 
et al. 2002), Maxima (Hanany et al. 2000), DASI (Halver- 
son et al. 2002), Cosmic Background Imager (Pearson et 
al. 2002) and Very Small Array (Scott et al. 2002) data. 
The 2dF galaxy redshift survey constraints are from Perci- 
val et al. (2002). The results obtained from the CMB data 
alone, using the 6 parameter model of Lewis & Bridle (2002) 
and fixing the optical depth to reionization r = 0.04 and 



Figure 10. The 68.3 and 95.4 per cent confidence contours on as 
as a function of n m from the present study (dark, solid curves; 
as in Fig. 3) together with the best-fit results from the cos- 
mic shear studies of Van Waerbeke et al. (2002: upper, dotted 
curve), Refregier et al. (2002: upper, short-dashed curve), Bacon 
et al. (2002: long-dashed curve), Hoekstra et al. (2002: dot-dashed 
curve), Hamana et al. (2002: lower, dotted curve), Brown et al. 
(2002: lower, dashed curve) and Jarvis et al. (2002: solid curve), 
A flat ACDM cosmology is assumed. 
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Figure 11. The 68.3 and 95.4 per cent confidence contours on 
as and n m from the analysis of CMB (outer, filled contours) and 
CMB+2dF data (solid lines). An optical depth to reionization, 
t = 0.04 is assumed. We have marginalized over all other param- 
eters (h, Q b h 2 and n). The inner, filled contours show the results 
from the present study (as in Fig. 7), including Gaussian priors 
of h = 0.72 ± 0.08, b = 0.93 ± 0.05 and n h h 2 = 0.0205 ± 0.0018. 
A flat ACDM cosmology is assumed in all cases. 
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marginalizing over h, Qbh 2 and n, are shown as the outer, 
shaded contours. The results obtained from the combined 
CMB+2dF data set, marginalizing over the same parame- 
ters, are shown as solid lines. A flat ACDM cosmology is 
assumed throughout. We see that the results from all three 
data sets are consistent at the 68 per cent confidence level. 
Our results are also in good agreement with the independent 
analyses of CMB+2dF data reported by Lahav et al. (2002), 
Percival et al. (2002) and Melchiorri & Silk (2002), the re- 
sults from CMB data and Type la supernovae presented by 
Jaffe et al. (2001), and the analysis of the IRAS Point Source 
Catalogue Redshift Survey presented by Plionis & Basilakos 
(2001). 



6 IMPLICATIONS FOR OTHER WORK 

The constraints on cosmological parameters reported here 
((T 8 = 0.695 ± 0.042 and fi m = 0.287 ± 0.036 for an as- 
sumed flat ACDM cosmology) are consistent with, though 
(in most cases) tighter than, those obtained from a number 
of other, recent studies based on the observed X-ray temper- 
ature and luminosity functions of galaxy clusters. (Schuecker 
et al. 2003 report results with similar precision to those pre- 
sented here). Our results are also consistent with current 
findings from studies of anisotropies in the CMB, the dis- 
tribution of galaxies in the 2dF galaxy redshift survey, the 
properties of type la supernovae, early results on large scale 
structure from the SDSS, and the most recent results from 
studies of cosmic shear. 

Our results have a number of implications for other cos- 
mological work. Firstly, the agreement between the cluster, 
CMB and 2dF results in Fig. 11 suggests that the optical 
depth to reionization is not large (r < 0.2; this issue will 
be examined in more detail in future work). Secondly, our 
result on as implies that the possible excess power detected 
at high multipoles (I ~ 2000-3000) in the CMB anisotropy 
power spectrum with the Cosmic Background Imager (Ma- 
son et al. 2002; Bond et al. 2002) is unlikely to be due to the 
Sunyaev-Zeldovich (SZ) effect. Our results also have impor- 
tant implications for future X-ray and SZ clusters surveys, 
since the number of clusters detected at high redshifts with 
high X-ray luminosities, temperatures and SZ fluxes will be 
much lower in a cr 8 ~ 0.7 universe than predicted by simu- 
lations with (78 ~ 1- 

Finally, we suggest that given the precision of the con- 
straints on Qm available from current Chandra / gas (z) data 
for relaxed clusters, and the relatively small systematic un- 
certainties involved, the fgas(z) data should be considered 
as a powerful probe in future cosmological work, comple- 
mentary to measurements of redshift evolution in the X-ray 
temperature and luminosity functions of galaxy clusters. 
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